Stability of fermionic Feshbach molecules in a Bose- Fermi mixture 
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In the wake of successful experiments in Fermi condensates, experimental attention is broadening 
to study resonant interactions in degenerate Bose-Fermi mixtures. Here we consider the properties 
and stability of the fermionic molecules that can be created in such a mixture near a Feshbach 
resonance (FR). To do this, we consider the two-body scattering matrix in the many-body envi- 
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ronment, and assess its complex poles. The stability properties of these molecules strongly depend 
on their centre-of-mass motion, because they must satisfy Fermi statistics. At low centre-of-mass 
momenta the molecules are more stable than in the absence of the environment (due to Pauli- 



somewhat less stable. 



blocking effects), while at high centre-of-mass momenta nontrivial many body effects render them 
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I. INTRODUCTION 

Nowadays the physics of cold atomic gases and their mixtures provide extraordinary 
opportunities both to test theoretical models and to predict new phenomena. The most re- 
markable effects arise when a resonant interaction can be engineered between pairs of atoms, 
by means of a magnetic-field Feshbach resonance (FR). The effects of resonant interactions 
are by now widely studied both theoretically and experimentally for both Bose and Fermi 
systems. The resonant interaction for two-component Fermi systems enabled the realization 
of BCS-BEC crossover (jj, y, La, |^ . The Bose system with a resonant interaction revealed 
such phenomena as the "Bose- nova" [5j]. 

To combine particles of both symmetries, and to consider a resonant interaction in a 
Bose-Fermi (BF) mixture, is a logical next step in this field. Recently, experiments have 
observed interspecies FR's between bosonic atoms and fermionic atoms p, LD, |8|, |9( . The BF 
FR brings new pos sibilities to the physics of q uantum gases such as boson-mediated Cooper 
pairing (13, [lj, lift , phase separation [ly, ll2| and the creation of a superfluid of fermionic 
polar molecules 18| . 

For trapped BF mixtures with a tunable interspecies interaction, one can expect a variety 
of regimes that are defined both by the interaction strength and by the densities of the two 
species. For example, the BF mixture may generate a new kind of a crossover upon crossing 
the resonance region. On the atomic side of the resonance, BCS superfluidity may be induced 
by the exchange of a boson density fluctuation, whereas superfluidity will likely disappear 
on the molecular side, as the molecules fill a Fermi sea. Further, if the crossing of the 
resonance can be accompanied by a process that produces ground state fermonic molecules 
(say, by a stimulated Raman process) then the resulting molecules would possess strong 
dipolar interactions that may restore superfluidity. 

But beyond this, the behaviour of the atom pairs in the resonant BF mixture depends 
subtly on the momentum of the pairs. Roughly speaking, for a single, free pair of atoms, 
a true molecular bound state exists on one side of the resonance, denoted as the "negative 
detuning" side. On the other, "positive detuning," side, the pair is not rigorously bound 
but may exhibit resonant scattering. For free molecules this demarcation at zero detuning 
between bound and resonant states is clearly independent of the pair's centre-of-mass mo- 
mentum. In the many-body environment, however, this situation changes dramatically. We 



will show that slowly moving molecules can be stable against decay even on the positive 
detuning side of resonance. The reverse is also true: pairs that are moving fast enough will 
become unstable and exhibit only resonances on the negative detuning side, even though 
their two-body analog would be completely stable. This unusual behaviour is connected to 
the fact that the pairs are themselves fermions, and must obey the correct Fermi statis- 
tics. This is of course different from the case of either a boson-boson or fermion-fermion 
resonance, where the pairs are always bosons. 

Both a quantitative and qualitative explanation of the resulting phenomena requires 
thorough implementation the two-body and many-body physics. Fortunately, a mean-field 
approximation is expected to give reasonable physical insight of physics near a FR JIOL IllL 



121 ] . In this article we will consider only one aspect of BF mixture near FR, namely how the 
stability of a composite fermionic molecule will be affected by the many-body medium. To 
do this we will assess the poles of the many-body T-matrix of a molecular pair propagating 
in the many-body medium. 

II. POLES OF THE T-MATRIX 
A. Two-body case 

The model we will use in the follo wing is somewhat well studied, having been introduced 



in similar contexts |26l . |27| : in Ref. |30j it was applied to the system at hand. We will 
therefore abstain from providing many of the details, and lean substantially on the results 
obtained in those works. 

The FR in two-body scattering is identified by parameterizing the dependence of the 
scattering length on the detuning, u, which represents the molecular binding energy in the 
case of an infinitely narrow resonance, and, in practice, is tuned using magnetic fields. We 
will then use the following parameterization for the scattering length 

a = aba ~ "2^-' (1) 

where ab g is the background scattering length away from resonance, m^f is the boson-fermion 
reduced mass, and g is a parameter representing the width of the resonance. For concrete- 



ness, we consider parameters appropriate FR that has been observed in in K-Rb 
Using the parameters defined in Ref. [30], we have a& 9 = —202 a , 5b = 5.1 x 10~ 5 K/G, 



and Ab = 1G, Bq = 513G, whereby the working parameters become v = 5b{B — Bq) and 



9 = JVbg8BA.B, where V bg = ~^ 3 - • In terms of these parameters, the low-energy two-body 



scattering T-matrix near a resonance takes the familiar form 
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m , ^ N 2nh 1 . , 

T 2B (E) = ^=, (2) 

m bf -a^ 1 + r m b fE - iJ2m b fE 



where ro = — 27r/m 2 jg 2 is the effective range of the interaction resulting from the underlying 
model. 

Bound states and resonances of the two-body system are identified in the structure of 
poles of (J2J in the complex plane. This is illustrated in Figure ^ where real and imaginary 
parts of the poles' energies are plotted as a function of detuning v. Here v is given in units 
of [j,, which is the same as the chemical potential in the many-body example below. We 
normalize v in this way for easier comparison between the two-body and many-body results. 
For v < 0, the two-body system possess a true bound state, whose binding energy is denoted 
by the solid line. In this case, the pole of (J2J) occurs for real energies. This bound state 
vanishes as the detuning goes to zero, which in fact is what defines the zero of detuning. 

For positive detunings v > 0, there is no longer a true bound state, but there may be 
a scattering resonance. This resonance appears for detunings v/li > g 4 mL-/47r 2 , about 
28 in Fig. ^ and its position is shown as a solid line, representing the real part of the 
pole. At the same detuning, there emerges a positive imaginary part of the pole energy 
(dashed line), which denotes the energy width of the resonance. Interestingly, for detunings 
< vj [i < g 4 ml r/4n 2 , the poles of (J2J are purely imaginary, and the imaginary part is 
negative. These poles stand for physically meaningless solutions to the Schrodinger equation, 
in which the amplitude in the resonant state would grow exponentially in time, rather than 
decay. These poles do not therefore identify any particular features in the energy- dependent 
cross section of the atoms. There is a characteristic detuning scale on which these events 
occur. This scale is given by the width parameter 7 2 = g 2 M b j /\/27r 23J] , also indicated in 
the figure. 

B. Many- body case 

In a many-body environment, the T-matrix and its poles depend on the centre-of-mass 
momentum. The importance of taking into account this dependence was demonstrated 
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FIG. 1: Complex poles of the two-body T-matrix, as a function of detuning. Solid and dashed 
lines denote the real and imaginary parts of physically relevant poles, respectively. The dotted and 
dash-dotted lines are real and imaginary parts of unphysical poles, respectively. 

in J20J for the BCS-BEC crossover at T ^ when non-condensed fermion pairs lead to 
pseudogap effects above T c and non BCS- behaviour below. The influence of the many- 
body medium on a Feshbach resonance between two fermions was demonstrated in [2l| for 
a composite boson as well. In the BF mixture the composite object is a fermion, so its 
momentum dependence cannot be omitted even for T = 0. Thus the T-matrix near the FR 
in a BF medium must be considered as a function of densities of both the Bose and Fermi 
subsystems as well as the centre-of-mass motion of a given BF pair. We have found that 
T-matrix poles and residues are quite sensitive to all these ingredients. 

The T-matrix T of the system in the many-body medium is defined by the Lippmann- 
Shwinger equation: 

f = gbg\ + gDg\G B G F f (3) 

where D is the renormalized molecular Green function, and Gb/f is the boson/fermion 



renormalized Green function, all of which are defined in terms of their well known 
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unrenormalized counterparts D°, and G°b/f by the self-consistent set of equations 

D = D° + D°flD 

Gb = G% + G%gG F Dg^G B 

G F = G° F + G F gG B Dg^G F 

where II is the molecular self energy, and the Green functions for fermions Gp, bosons Gb 
and molecules D represent 2 ® 2 functions containing their anomalous (pairing) parts. All 
quantities are functions of the energy E and centre-of-mass momentum P of the molecules. 
The complete solution to these equations is beyond current computational capabilities, ex- 
cept perhaps by Monte Carlo methods. 

We therefore make a few simplifying assumptions, namely we account for propagation of 
the atomic fermions and bosons using their free Green functions only (i.e., setting Gp ~ G F 
and Gb ~ G B )• An important consequence of these choices, is that the many-body T- 
matrix is approximated by its ladder series, which means that II m gG B G F g^. This standard 
approximation has the property of being exact in the 2 body limit, where Eq. (J2J) satisfies 
(JHJ). This implies that the two body physics is accounted for exactly in the many body 
problem. 

Another appealing characteristic of this approach is that the molecular self energy II can 
be calculated exactly J2J|, leading to the following expression 

U(E, P) = £m bf A - £m bf k f - & (^ 

_ m lf P _ m^A i ( (k f +P m bf /m b )*-D \ 
m b P J lLL \(k f -P m bf /m b y-D) 

+ 4^ m bfV- U lZl \(k f -jDF-(P m bf /m b r J > W 



where D = J2rribf(E — P 2 /2(rrif + m h ) + //), and A is an ultraviolet cutoff, which can be 

2 

regularized by a shift in the detuning v — > v — ^rribfA. 

Further approximations include considering a homogeneous system where the density of 
the fermionic subsystem is much smaller than that of the bosonic one. This is important 
for treating the bosonic subsystem as almost unperturbed by the fermionic one. Finally, we 
disregard the boson-boson interaction, and include only the resonant part of the Bose- Fermi 
interaction, setting a bg = 0, in Eq. JTJ), though not in the definition of g. 



The set of equations we obtain with these approximations is therefore 

T(E)=g 2 D(E), (5) 

where D(E) is the pair propagator defined by: 

D(E)=D°(E) + D°(E)U(E)D(E), (6) 

U(E)=g 2 G° B (E)G F (E), 



which lead to 



9 2 



T(E, P) = - 2 2 . (7) 

In order to understand the stability of a BF molecule, we study the structures of the 
poles and residues of the T-matrix in Eq. (j7j) . To determine numerically the poles, we 
consider this equation as a system of two nonlinear equations for real and imaginary parts 
of the internal energy, E — P 2 /2(rrif + m^) + \x. The resulting nonlinear system of equations 
always has some unphysical solutions which can be rejected by the following analysis: i) the 
residues for physical solutions must be less than unity; ii) for the imaginary solutions the 



relative momentum, D = J2rribf(E — P 2 /2(m/ + m&) + //) in Eq. PJ), should dwell on the 
lower half of the complex momentum plane (the second sheet for a complex energy), as for 
resonant scattering [28(; iii) the sum rule (fTTl see below) should be fulfilled including both 
the discrete and continuum parts. 

As an example, we have calculated the poles of (J7J) for a dual species gas with fermion 
density 10 13 cm -3 and boson density 10 14 cm -3 . Moreover, for simplicity, we assume that 
the bosons and fermions have the same mass, m = 40 amu. The results are shown in Figure 
El for two representative centre-of-mass momenta P of the molecules. (Note the different 
scale from Figure HJ. Unphysical poles are not shown in this figure. These results are cast 
as a kind of "binding energy," by subtracting the centre-of-mass kinetic energy, and adding 
the chemical potential \i. They can therefore be compared directly to the two-body results 
in Fig. [D 

The structure of these poles is quite different near resonance, although we stress that far 
from resonance [y ^> 7 2 /2) they return to the two-body values. A main difference from 
the two-body case is that now there may exist physical poles for any detuning across the 
resonance. The appearance of poles depends, however, on the centre-of-mass momentum P 
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FIG. 2: Complex poles of the many-body T-matrix, as a function of detuning. Solid and dashed 
lines denote the real and imaginary parts of physically relevant poles. In (a) the centre-of-mass 
momentum of the molecule is P = kf, i.e., equal to the Fermi momentum of the atomic gas. In 
(b), this momentum is P = O.lkf. 



of two atoms. If P ^> 2kf, the two-body physics will not be influenced much and the poles 
structure will be as it is in Fig. ^ We have verified this numerically, but have not shown it 
in this paper. 

Figure I21^a) shows the case where the centre-of-mass momentum of the atom pair is equal 
to the Fermi momentum of the atomic gas, P = kf. The criterion of molecular stability is set 
by the vanishing of the imaginary part of the pole. Figure|2ta) shows that the imaginary part 
remains zero until u/fi ~ —5, after which the pair becomes unstable. Thus, for detunings 
—5 < u/fi < 0, some of the molecules that would have been stable become de-stabilized 
in the many-body environment. They may still be rather long-lived, however. In the limit 
of very large momenta, P ^> 2k f, we find that the pole structure returns to the two-body 
value. This makes sense, since a very rapidly-moving molecule does not interact strongly 
with the gas at all. In particular, for P > 2k f, we find that a "gap" reappears, in which 
no physical poles exist for some range of positive detuning. In the P — ► oo limit, this gap 
returns to its two-body value, 7 2 /2. 

Figure Efb) shows an alternative case in which the molecular momentum is much smaller 
than the atomic Fermi momentum, P = O.lkf. In this case, the imaginary part of the 
pole only differs from zero at positive detunings, vj ' \i > 8. Therefore, quite the opposite 
to the P = kf case, here the molecules that would have been unstable are stabilized by 
the presence of the many-body environment. Roughly, this is due to Pauli blocking of the 
fermionic atoms into which the molecules would dissociate. If such an atom already occupies 
the state into which the molecule would drop its atom, then the process is forbidden. We 
discuss this further below. For any molecule with P < 2k f, there exist physical poles at all 
detunings. Thus positive-detuning molecules are always present in the BF mixture. 

III. CONDITIONS FOR MOLECULAR STABILITY 

To make a more global picture of whether molecules are stable or not for a given detun- 
ing, we can consider the spectral function for atom pairs. In general, when the molecular 
propagator possess an imaginary part, this part alludes to the decay rate of the pair due 
to interactions with the rest of the gas, in the same sense that the oscillator strength of an 
atom alludes to its decay rate by spontaneous emission. Following a standard approach [32] , 



the spectral function is defined as 

p(E, P) = --ImD(P, E). (8) 

7T 

In the case of a true, bound molecular state, the spectral function reduces to a delta function 
at the energy of the state E : 

p(E,P) = 2vr Z(P)5(E-E ), (9) 

where the coefficient is given by the "spectral weight" function 

Z(P) = l-dRe(Il(P,E))/dE\ Eo ' (10) 

Just as for an oscillator strength, the discrete and continuum parts of the spectral density 
must satisfy a sum rule: 

Z(P)+ JdEp(E,P) = l (11) 

for each momentum P. We have explicitly verified the sum rule in each case we computed, 
as a test of the numerical procedure, and to distinguish between physical and unphysical 
poles of the T-matrix. 

As mentioned above, the spectral weight Z is associated with a specific pole E of the 
T-matrix, an energy level of the system, and represents its population. Therefore, if Z 
vanishes, so does the probability of finding stable molecules in the gas. To this end, Figure 
IIHI plots contours of the function Z(P) as a function of detuning v and centre-of-mass 
momentum P of the atom pairs. These calculations are performed for fermion and boson 
densities nf = 10 13 cm -3 and n^ = 10 14 cm -3 . The contour of Z = thus represents the 
borderline between conditions where molecules exist and are stable (to the left of this line) 
and where they are unstable to decay (to the right of this line, in the white region of the 
graph) . 

Figure IIHI thus shows that: i) molecules are still stable for a continuum of positive 
detuning when P is small; ii) molecules that would have been stable at negative detuning 
may not be such stable at intermediate momenta P (although at small negative detuning 
they may possess small widths); and iii) in the limit P — > oo, the borderline between stable 
and unstable again returns to zero detuning. Detailed information on the spectral density 



can also be used to elicit t 
completed in future work 



re momentum distribution of molecular pairs, a task that will be 
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FIG. 3: Contours of molecular population under various combinations of molecular centre-of-mass 
momentum P and detuning v. The uppermost contour identifies the detuning at which bound 
molecules first appear for a given momentum P. Numbers indicate contours with a equal molecule 
fraction. The dots represent the result obtained analytically for the critical detuning, Eq. (fTB|) . 

Thus far, these are rigorous results, at least within the simplifying approximations made 
above. Once this is done, the effective dissociation energy of the molecules within the 
medium is determined. The relation between the molecule's total energy at dissociation 
and the molecule's momentum is then easily determined from kinematics, plus simple con- 
siderations on the Pauli blocking of the atomic fermions. For example, consider the case 
where the molecule's kinetic energy is greater than twice the atomic Fermi energy, i.e., 
P 2 /2(mf + nib) > 2 x k 2 J2m. At the same time, the molecule is assumed to be exactly 
at its dissociation threshold, so that it could live equally well as a molecule or as two free 
atoms. Upon dissociating, each atom would carry away half the energy, so in particular 
the fermionic atom is at the top of the Fermi sea, and this dissociation is not prevented by 
the Pauli exclusion principle. The total energy of the molecule at its dissociation threshold 
is then determined simply by the molecular kinetic energy, and no contribution is required 
from the molecular binding energy. Thus, if B represents the internal energy of the pairs 
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relative to threshold, bound molecules are possible when 

B < for P > 



2(mb + m V (12) 



Alternatively, suppose the molecules have less than twice the atomic Fermi energy, 
P 2 /2(rrif + rrib) < 2 x k"j/2m. Now it is no longer guaranteed that the molecules can 
automatically decay in the many-body environment, since the fermion's kinetic energy may 
lie below the Fermi level of the atomic gas. In such a case, the molecule can sustain a 
positive internal energy without dissociating, simply due to Pauli blocking. To decide how 
high this binding energy can be, we examine the conservation of energy and momentum in 
the dissociation process: 

P = P/ + Ps (13) 

p2 . b = jL + jL. (14) 



2(2m) 2m,b 2m/ 

Here pj, p^, and P are the momenta of the atomic fermions, atomic bosons, and molecules, 
respectively. To ensure that the atomic fermion emerges with the maximum possible kinetic 
energy, we consider the case where P and p/ point in the same direction. To ensure that 
Pf > kf, where kf is the atomic Fermi momentum, along with ([14)1 . implies that molecules 
are stable when 



B < (m f P-(m b + m f )k f Y ^ p < 



2 ^ + m f\ kf . (15) 



2rrib vrif (rrib + mf) \ m/ 

Figure 0] shows the internal energy of the molecules evaluated at the stability boundary, 
as described above, as a function of centre of mass momentum. The solid line in this figure 
is determined numerically from the Z = contour of Fig. 11111 Subtracting the kinetic energy 
contribution and chemical potential from the pole of the T-matrix evaluated on the contour, 
we obtain the molecular internal energy. Also shown, as dots, are the kinematic estimates 

(gasp . 

Analytical expressions for the detuning at the boundary as a function of centre of mass 

momentum, plotted as dots in Fig. IIH[ are readily obtainable analytically, with similar 

accuracy by inspecting the denominator of Eq. (JJJ). The total energy of the molecule, 

measured from the chemical potential, is, in fact, given by the pole of ((7j). In general this 

p 2 

2(m,f+m b ) 



energy can be written as Eq = », p — r — /i + B, where B is a complicated function of all 
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FIG. 4: The internal energy that a molecule of centre-of-mass momentum P/kf would require to 
overcome Pauli blocking and dissociate. The shaded area represents numerical calculations, and 
the dots represent the analytical kinematic arguments in Eqs. (|12I15|) , 



the parameters. However, since E to t is a pole of (jZJ), then E Q 

so 

/ P 2 

Vcrit — B — IT 



-fi + B,p). (16) 



2{m,f + m b ) 

Plugging the stability boundary value of B from Eqs. ()12|15j) into this formula leads to an 
analytic, albeit complicated, expression for the critical detuning as a function of centre of 
mass momentum. 



IV. CONCLUSIONS 

The constitution and stability of a composite fermion in the many-body environment is 
an important part of the BF crossover regime. Within our approach we have addressed 
properties of strongly correlated BF pairs and defined the stability region. We concluded 
that in the many-body environment, low momentum molecules can be stabilized at positive 
detuning, and intermediate-momentum molecules can be de-stabilized, existing for shorter 

13 



times at small negative detunings. We also concluded that there is always a probability 
to observe molecules at positive detunings, provided their momenta are less than « 2pj, 
even though these molecules would not exist as two-body objects. The way in which these 
molecules are distributed will form the basis of future work. 
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